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Abstract. This is a series of lectures on Monte Carlo results on the non-perturbative, lattice 
formulation approach to quantum field theory. Emphasis is put on 4D scalar quantum field theory. 
I discuss real space renormalization group, fixed point properties and logarithmic corrections, par- 
tition function zeroes, the triviality bound on the Higgs mass, finite size effects, Goldstone bosons 
and chiral perturbation theory, and the determination of scattering phase shifts for some scalar 
models. 
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1. Lattice Quantum Field Theory 

In this series of lectures I attempt to cover some numerical applications to stochas- 
tic analysis in relativistic quantum field theory (QFT). Following a more general 
introduction I then focus on the non-perturbative quantization of scalar QFT. This 
emphasis is dictated by the limited time and place available as well as by the progress 
in the subject. Scalar QFT meanwhile is the most advanced of the interesting mod- 
els, both in theoretical understanding as well as in numerical stochastic calculations. 

If you are interested in particular topics of lattice field (including gauge field) 
theory, a good starting point is to look up the proceedings of the annual workshops 
in lattice field theory. For the recent years they have been published in Nucl. Phys. B 
(Proc. Suppl.); up to now there are volumes 9(1989), 17(1990), 20(1991), 26(1992), 
30(1993). In every one of these volumes there are also reviews, although mainly 
targeted at the devotees of the field. 

There are only few books giving surveys. Almost all important results of the first, 
hot days of Monte Carlo computations are given in the reprint collection [116] and 
the monograph [24]. The more rigorously oriented person should maybe start with 
[120]. Recently two volumes on the field have been published, one is the monograph 
[118] and the other contains a collection of review article written by experts on 
various subtopics [25]. 

In the first part I attempt an overview on lattice field theory and the related 
computational issues. The role of the model in particle physics is discussed and 
the lattice formulation for that model is introduced. In the second chapter some 
rigorous statements on the model are reviewed and then lattice calculation results 
concerning the fixed point structure and the issue of triviality are presented. In the 
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subsequent chapter further results on the so-called triviality bound are discussed; 
also the particle spectrum and in particular the finite size efi^ects produced by the 
massless Goldstone states are reviewed. In the last section I present recent results 
that use the scalar model as a testing ground for new methods to determine scattering 
phase shifts non-perturbatively. 



1.1. Regularization and Quantization 

1.1.1. Space-time Lattice 

Quantum field theory usually is discussed within the perturbative framework. Even 
there it is obvious that there arc divergences in elementary contributions. In a func- 
tional integration (Feynman path integral) formulation one finds that the integration 
is not well-defined. The functional integral has to be regularized in order to be able 
to identify the singularities. 
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Replacing the space-time continuum by a lattice provides a suitable regulariza- 
tion as well as a computational tool [142]. We can define the continuum QFT as 
the limit of the lattice QFT for vanishing lattice spacings. Important symmetries of 
the Lagrangian (like e.g. gauge invariance) may be sustained in the lattice formu- 
lation. On the other hand, obvious symmetries of space-time like translational and 
rotational invariance are replaced by their lattice subgroups and eventual recovery 
in the continuum limit has to be established. 

We will discuss euclidean QFT. The regularization replaces the continuous, 0(4) 
symmetric space-time arena by a (in practical calculations) finite lattice with sym- 
metry 0(4;Z). 
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lattice A e.g. c 

field $ : A X (e.g. C) G T 

configuration F e 

action S'a : F i-^ K 

partition function Za = / (i/iA(*i') exp[S'A((7, 'J')] 

free energy J^a = In ^A 

Gibbs factor PAi^) = ^a ^ exp[S'A(.g, $)] 

expectation values (^)a = / d^A{^)pAi^)A{^) 

Quantization thus amounts to determination of expectation values in a 4D statistical 
system, like e.g. the Ising model (cf. for instance [117]). With g we denote the 
coupling constant (or a vector of constants) parametrizing the action. 

1.1.2. Scaling 

The focus of interest in that process are two limits: 
Thermodynamic limit: The infinite volume limit |A| — > oo. 

Continuum limit: Is obtained at a critical point g gc in coupling constant 
space, where the lattice spacing constant a{g) — s- 0; the lattice spacing acts as 
a regularization parameter. 

The second limit needs some discussion. The physical quantities like mass or length 
of the continuum theory should approach constant values, however on the lattice 
they are expressed in multiples of the lattice spacing a. These lattice quantities are 
dimensionless and are determined from the connected correlation functions (propa- 
gators) of some operators 

{A{0)A{x))c oc exp(— a;/^) = exp(— ax/a^) 

= exp(-mphysaa;) , i.e. Wphysfl = 1/C • (1-1) 

These correlation functions describe propagation in euclidean time and the operators 
may be representing n-particle states. The exponential decay should hold for sensible 
theories [120]. „ 

Only at a critical point the correlation length divergesEl and allows to let a ^ 
while keeping a fixed value for the physical mass parameter. One has to fix the scale 



^ At the critical point the correlation function decay polynomially oc ^+'' 
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a{g) by trading one physical parameter. Once the scale is set, all other parameters 
may then be expressed in terms of this scale parameter a{g). However, all other 
physical parameters should scale correspondingly, only then can we identify the 
correct continuum theory. 

In practical situations we distinguish between two qualities of scaling behaviour. 

Scaling: The dimensionless ratios like — scale towards the physical 

ratio — as a ^ 

m2 

Asymptotic scaling: The scaling behaviour a{g) cx /(<?) as determined from the 
renormalization group /3-function. 

The relative importance of irrelevant operators determines the quality of scaling. 
Whereas scaling may be observed quite early one has to be quite close to the critical 
point in order to observe asymptotic scaling. 



1.1.3. Continuum Limit 

Let us collect some points important in the context of the continuum limit. 

— The lattice action should obey Osterwalder-Seiler positivity when approach- 
ing the continuum limit in order to guarantee existence of a transfer matrix 
formalism and positive norm Hilbcrt states [115]. 

— The continuum limit is obtained at a critical point (2nd order phase transition) 
of the lattice system: g ^ g* , where the dimensionless correlation length ^ — > 
oo. 

— All physical quantities are multiples of the lattice spacing and should scale 
correspondingly, a{g) — > for — > gi*. 

— Lattice translation invariance mod(a) should become continuous translation in- 
var iance. 

— The space-time symmetry (cubic group for cubic lattices) should become the 
continuous 0(4) (rotational) invariance. This may be observed in off-shell ob- 
servables (like the 1-particle propagator) as well as in on-shell ones (like differ- 
ential scattering cross-sections). 

— Fields, couplings and masses in the lattice action are bare quantities and are 
renormalized in the quantization in continuum limit. Renormalized quantities 
will have a subscript r henceforth. 

— Because of renormalization, the continuum properties of a quantized theory 
may be quite different from those of the bare theory. In statistical physics 
this observation is linked to the notion of universality: Theories with different 
micro-structure may show universal critical (macroscopic) properties. Since the 
critical properties are just the physical properties of our quantized continuum 
theory this limits the class of liable continuum QFTs. 

All these points should be observed in a sensible investigation of the continuum 
limit. Most of these properties can be studied only in the quantized theory in a 
non-perturbative environment. In existing Monte Carlo calculations several of these 
features have been checked and some of these checks will be discussed in this series 
of lectures. 
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1.2. Particle Physics and the Model 

1.2.1. Lattice Fields 

Typically we have to deal with the following types of lattice fields. 

Scalar bosons: Occur in the action usually with local terms like ^^'^ 
nearest neighbour terms like ^xUx.fj.^x+ti', the n.n. term is a lattice discretiza- 
tion of the continuum derivative term. Monte Carlo integration methods for 
these fields include the standard Metropolis or heat bath updating as well as 
cluster- and multigrid methods. The computational effort grows with 0(i*) up 
to 0{L^) depending on the algorithm {L denotes the linear extension of the 4D 
lattice), 

Gauge (vector) bosons: Action terms include the Wilson plaquettc interaction 
TrUp, where the plaquette variable is the path ordered product of link fields 
Ux^ij, around a plaquette. The gauge field coupling /3 cx where g is the 
usual perturbative coupling constant. Monte Carlo updating algorithms include 
Metropolis, heat bath and overrelaxation; some progress has been made in 
designing multigrid algorithms. The computational effort grows with 0{L^) 
up to 0{L^). 

Fermions: Occur in the action with terms like ^x'^x and ^xQxy^y which includes 
Dirac 7-matrices and nearest neighbour coupling k. Updating algorithms are 
an order of magnitude more demanding with regard to computational resources; 
they include molecular dynamics updating, Langevin and "hybrid" Monte Carlo 
methods as well as some recent attempts towards multigrid methods. The 
central numerical problem is the inversion of the lattice Dirac operator Q, thus 
the total effort grows with 0{L^) up to 0{L^^) depending on the algorithm. 

For the formulation of the models on finite volumes one usually assumes periodic 
boundary conditions, <&(.t + Lfi) = $(a;) (for fermions: antiperiodic b.c). The 
implicit assumption is that the thermodynamic limit should be independent on the 
choice of boundary conditions. Periodic b.c. seem to produce relatively small finite 
size effects. 

1.2.2. Scalar Model 

The standard model of electroweak and strong interactions describes extremely suc- 
cessfully nowadays results in particles physics in an energy domain from a few MeV 
up to more than 100 GeV. Its central symmetry is 5'C/(2) x U{1) x SU{3). The 
SU{3) is that of QCD, represented by quarks, gluons and their interaction. The 
non-perturbative lattice study of that model is a central topic in lattice calculaticms 
[29], but I will not try to cover it here. The symmetry groups SU{2) x [/(I) arc 
those of quantum fiavour dynamics and are represented by the fundamental Icptons, 
4 gauge bosons and the dublet $ = ($1 + »$2, ^3 + «$4) of complex Higgs fields. In 
particular the action for the Higgs fields has a symmetry SU{2)l x SU{2)ji isomor- 
phic to 0(4). Thus we may represent the Higgs fields by a vector of 4 real scalar 
fields, too. 
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Fig. 1. Standard model symmetry breaking scenario. 



In the standard scenario one assumes that there is a symmetry breaking potential 
hke the well-known Mexican hat potential A(<I>^ ~ 1)^- In the Goldstone realization 
the vacuum chooses spontaneously one direction (out of the infinite set of inequiv- 
alent representations, cf. fig. |^), thus one is left with one massive component (the 
Higgs particle) and three massless components (the Goldstone modes), which, in the 
so-called Higgs mechanism, become the longitudinal modes of three of the (originally 
massless) vector bosons. The three vector bosons acquire mass in that process. For 
sake of the later discussion let us distinguish the two processes: 

Symmetry breaking: The 0(4) symmetry is spontaneously broken and in the new 
vacuum we find one massive and three massless modes. 

Mass generation: Three massless gauge bosons buddy with the three Goldstone 
modes and become massive gauge bosons, Z and . 

Both mechanisms are non-perturbative as they involve infinitely many degrees 
of freedom. The results mentioned are based on tree-level perturbation theory and 
some additional assumption about the vacuum structure. One of the basic questions 
to be answered in the non-perturbative lattice simulations is, whether these results 
survive quantization and whether they remain valid at arbitrary Higgs coupling 
strength. 

The essential non-perturbative ingredient is the 0(4) model. This is the main 
motivation for studying that theory non-perturbatively in some detail. 

1.2.3. The ^'^ Model 
The continuum action 



has three parameters, the bare mass, the 4-point coupling, and the external field. 
Discretization is not unique and a particular simple form of a lattice action is ob- 
tained by replacing the derivative term by a nearest neighbour difference. With 




(1.2) 
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suitable redefinitions of the couplings we get 

Slatt = [-2ti J2 + - 1)' + J^x] ■ (1.3) 

A 

We have several important cases. 

— <i>a; S Z2: the Ising model, 

— <i>a; S K: the simple 1-component ^'^ model, 

— = €: 0(4): the 4-component model, i.e. the Higgs model. 
For notational simplicity we sometimes use field variables 

(p EE \/2^$ , j = J/V2^ . (1.4) 



In the form (L3) we identify immediately two limits of the system. For A — > we 
approach a gaussian model. For A ^ 00 only configurations with — 1 contribute 
and we have a nonlinear sigma model (for the 1-component version: Ising model). 
The coupling k, A and j in the lattice action are bare quantities, given at the generic 
lattice scale a. Values measured in the quantized system like e.g. the mass or the 
4-point coupling are renormalized quantities and will be denoted by a subscript: 

Ulr , At- . 

At the critical point of the model, k = Kc, critical exponents characterize the 
singular behaviour of various quantities [125, 104]. Introducing the reduced coupling 

t^^^ (1.5) 

we summarize in table | the critical behaviour and definitions of the critical indices. 
In sect. 2 we discuss the issue of triviality of the model. If it is indeed trivial (as 
all evidence points to) the critical exponents are those of mean field theory and 
there are non-leading logarithmic corrections. The exponents of these corrections 
are given in table | and depend on N of the group 0{N) [14, 83] {N = 1 for the 
one-component model). 

The two definitions given for the renormalized charge agree in the continuum 
limit (sometimes one uses the continuum definition action with an extra factor 6: 
\r = 3(mr / (fr))'^) ■ Here Z denotes the field renormalization constant (and will be 
discussed in sect. 3 in more detail) and not the partition function (sorry, but thats 
how it is). Only two of the critical exponents are independent (cf. e.g. [76]). The 
mean field values are for the scaling exponent of the correlation length — ^ and for 
the anomalous dimension rj = 0. From the scaling relations we find a — 2 — di' — 0, 
S = ^ = 3, 7 = K2 - ^) = 1 and /? = j/{S - 1) = i 

1.3. Monte Carlo Methods 

Quantization amounts to determination of expectation values 

(^)-|E^We^p[^(*)]- (1-6) 
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TABLE I 

Summary of the expected critical behaviour in the AT-component model in terms 
of the reduced critical coupling t = (kc — k)/kc and the external field j. 



Order parameter 

(*> ~ Itl^llnltll"^ forj = 



yl/^llnljlli fort = 



Susceptibility 

X=($(p = 0)2>e = ^^(*0*x>c ~ |t|-T|ln|t||wil 



Internal (link) energy 

{E) = ^{*,*,+^> 



Specific heat 

= '"'^ = JX\^^^'^ ' ^^^'^ ^ ir"iin|tii^ 

Correlation length 

(^o^x)c oc exp(-TOr|x|) for t ^ 

oc l/\x\'^-^+'^ for t = 
mr = l/(, ~ Itl^llnltll^W+ie 



Renormalized 4-pt coupling 

(^(p = 0)4), 

^{*(p = 0)2)2 

or = -imr/i'Pr)? ~ llnKir^ 

where (i^r-) = iv)/^^ 



-N-2 



In Monte Carlo integration one samples the field configurations with an equilibrium 
probability distribution -P($) — exp[S'($)]. This allows to determine expectation 
values according 

{A) = lim 1 / rfr A($(r)) = lim 1 A($„). (1.7) 

T— >(x> 1 Ici iv— >oo iV '—^ 

n=l 

The stochastic variable is the field $(t) evolving with the stochastic (computer) time 
r. It is determined from a homogeneous Markov process with transition probability 
T, 

P{^n = Vn\^n-1 = ^n-l) = Ti((^„|<p„_i) (1.8) 
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i.e. Tfc depends only on the numbers of time steps k but not on the history. If 
Tk{ip'\(f) > for all pairs {(p',f) and the process is aperiodic then it is ergodic. 
A Monte Carlo integration algorithm usually constructs such a Markov sequence of 
configurations, similar to a sequence of snap shot of some realistic ferromagnet. Thus 
the reason for calling the method a simulation. In realistic simulations ergodicity 
and the related problem of relaxation are important questions. In particular if there 
are topologically different sectors in configuration space some Monte Carlo updating 
algorithms may have problems connecting them. 

If implemented correctly the Monte Carlo integration leads to correct results with 
a purely statistical error that decreases with the number of iterations like 1/\/N. 
The results are for finite size lattices, however. Control of both, the statistical error 
and the finite size effects, belongs to the central issues of this field. 

One starts with an (arbitrary) configuration. The Markov sequence of configu- 
rations 

(po ^ ^ ^2 ^ ■ ■ ■ (1-9) 

produces samples of the corresponding probability distributions, which, depending 
on certain sufficient conditions, approach the equilibrium distribution, 

p(0) ^ p(i) ^ p(2) _^ (equilibrium distribution), (1.10) 

where P(") = r"p("). In a calculation one therefore starts to measure observables 
only after a sufficient number of "equilibrating" Monte Carlo steps. Subsequent 
updates then produce configurations according to the required distribution. 

Meanwhile wo have a collection of updating algorithms [145]. Some of them 
proceed by changing the configurations locally, i.e. just one variable at a time, 
others try to implement collective updates of nonlocal nature. 

Local updating algorithms 

— Metropolis: The simplest and oldest [107] algorithm: we discuss it below. 

— Heat bath: Like an iterated Metropolis updating, optimizing the local 
acceptance rate. 

— Overrelaxation: A sometimes very effective generalization of the heat 
bath algorithm. 

— Micro canonical: Reformulates the updating in terms of a deterministic, 

discrete Hamiltonian evolution in a space with doubled number of vari- 
ables: ($,n). A variant method introduces a microcanonical demon. 

— Hybrid Monte Carlo: Combines a number of microcanonical updates 
(a "trajectory") with a final Monte Carlo acceptance step. 

— Langevin: Uses the stochastic differential equation for the construction of 
configurations; equivalence to the microcanonical and the hybrid method 
may be demonstrated (in certain limits). 

— Multicanonical: Or multimagnetical or entropy sampling, samples an 
ensemble with modified energy density, particularly well suited for systems 
with degenerate ground states (like at 1st order phase transitions or in 
spin-glasses) [5, 6]. 
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Non-local updating algorithms 

— Multigrid: Employs scale changing transformations to identify possible 
collective updating modes, essentially a simultaneous updating on several 
length scales. 

— Cluster: Identifies the system-specific dynamical nonlocal degrees of free- 
dom; this method is very efficient for certain models and will be discussed 
below in more detail. 



1.3.1. Metropolis Algorithm 

Considering a stochastic process X{t), detailed balance is a sufficient condition that 
the probability distribution approaches the equilibrium distribution exp [5'(a;)], 

Pi(x|2/)e^(^) = Pi(y|x)e^(") ^ lim P(X(r) = x) ^ e^(^). (1.11) 

r— >oo 

Summation over x shows that the limiting distribution is an eigenvector of Pi . The 
Metropolis algorithm [107], starting at some value x (i.e. some initial configuration), 
consists of the steps: 

1. Choose some y according to some a priori transition probability PT{y\x) . 

2. Accept the new value y with the probability 

PA{y\x) = min ( 1, ^^N,) ^^Pjf (^jh (1.12) 
= min(l, exp[S(j/) - S{x)]) for PT{x\y) = PT{y\x) . (1.13) 

It is straightforward to see, that the resulting Pi = PtPa fulfills detailed balance. 
In particular for symmetric Pt the information necessary to decide on acceptance 
or rejection comes only from the change of the action AS" with regard to the 
change of the variable. If this change is local, i.e. just at some field variable 
the change of action may be determined from the field values in the local 
neighbourhood. 

3. Repeat these steps from the begin. 

One MC iteration is one sweep through the lattice, that is one attempt to update for 
each field variable on the lattice. Typical calculations perform O(IO^) MC iterations 
for given lattice size and couplings constant values. 

The exponential autocorrelation time controls the statistical dependence of sub- 
sequent configurations and may be defined 

|p(n+l) _ p| < exp (-l/Te,p)|P(") - P| (1.14) 

or from the autocorrelation of some observable A, 

VA{n) = {A{ipQ)A{iPn)) - {A)"^ ~ exp {-n/Texp,A) , Texp,A = SUpTexp,A , (1-15) 

A 

or from the integrated autocorrelation time 

n„M = 2 + Z.r40) • (^-^^^ 

n>0 ' 
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The latter quantity controls directly the statistical error in estimates of (A); a run 
of n measurements contains only n/2Tint,A effectively independent measurements. 
This relaxation behaviour affects particularly systems with large correlation lengths 
5 (near critical points) , since the autocorrelation length grows like 

Tea;p ^ min(i,^)'= , Tint^A ^ min(i,^)'='^ (1-17) 

where z and za denote the dynamical critical exponents. They depend on the 
updating algorithm and have the value 2 for most local algorithms like the Metropolis 
algorithm (for overrelaxation algorithms this value may be smaller: ^ 1 [145]). The 
computational effort for the simulation of a system of linear extension L near a 
critical pointa, where ^ ~ L, grows like L'^'^^^. This phenomenon is called critical 
slowing down. For more detailed discussions cf. [123, 124, 145]. 

In the MC simulation the statistical error of the measured quantities shrinks 
with the inverse square root of the number of statistically independent configura- 
tions. Improvement of the performance therefore can be achieved only by improving 
the updating algorithm by making them computationally simple (and quick), by 
reducing the autocorrelation time and by reducing the dynamical critical exponent 
z. This is the motivation for the introduction of multigrid- and cluster updating 
algorithms. 

1.3.2. Cluster Algorithm 

Fortuin and Kastelyn [81, 82] pointed out that the Ising model has an alternative 
formulation as a bond percolation model. This leads to clusters of connected bonds 
[23], which are the nonlocal dynamical objects of the Ising model. 

Swendsen and Wang [131] suggested a very efficient updating algorithm based on 
alternating between the spin and the cluster representation. For a given spin config- 
uration one first identifies Ising clusters built of neighboured sites with parallel spins. 
Some of the bonds are then deleted with probability exp(— 2k), the resulting bond 
clusters are the dynamical FK-CK-SW-Clusters [81, 23, 131]. With a probability 0.5 
one now can flip simultaneously all spins in a cluster, repeating this for all clusters. 
This results in a new spin configuration and the updating step is completed (cf. fig. 
i- 

It turns out that this algorithm gives z ~ 0.3 for the 2D and z ~ 1 for the 4D 
Ising model. Wolff proposed a modification, where one randomly chooses a site and 
identifies the cluster connected to it like in the SW-algorithm; this cluster is then 
flipped unconditionally. This results in z ~ for the 4D system! 

The cluster algorithm originally was designed for systems like the Ising or the 
Potts model, where the spin variables take values in a discrete group. Wolff [144] 
observed, that one may extend its scope of application to continuous groups, that 
have embedded Ising variables (the group manifold may be split into disjoint parts, 
cf. [124]), like e.g. the group 0{N). This allowed to utilize the updating method for 
A^-component $^ theories. The cluster method has provided a major breakthrough, 
the efficiency of the algorithm has brought forward results for scalar models with 
highest computational precision. 

^ For systems with several coexisting phases the autocorrelation length grows even exponentially 
with ~ exp(cL'*~^); multicanonical algorithms address this tunneling problem. 
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Fig. 2. The Swendsen-Wang cluster updating algorithm: {s} {b} —* {b'} {s'}. 

The cluster representation allows the construction of observables with reduced 
variance, so-called improved estimators [127]. The propagator between spins at two 
sites, e.g., may be obtained by just counting how often both spins belong to the 
some cluster: 

= ( Xc{x)Xc{y)) where Xc{x) = | ^ ![ ^ ^ ^ . (1.18) 

all clusters c 

For a given set of configurations this allows determination of the propagators for 
larger distances with much better statistical quality than achieved from the direct 
measurement. 

The clusters are the dynamical objects of the critical system, a realization of 
Fisher's droplets. One therefore would expect relations between the geometrical 
properties and some critical indices. Due to scale invariance at a critical point 
clusters should look alike on different length scales. In quantum mechanics this is 
expressed by the fractal structure of the quantum paths, which are of Brownian 
nature. We might expect such a fractal structure also in the field configurations of 
quantum field theory (cf. fig. 

The relation between such a fractal structure and the odd critical index may be 
argued as follows [126]. Under a scale transformation we have scaling of the reduced 
critical temperature t (cf. table ^ in the form 

t^h-^'^t (1.19) 

where v denotes the scaling index of the correlation length. In a volume L'^ the total 
magnetization scales according to table | 

M{L) cx LHI^ -> (1.20) 
M{hL) cx [hLYh-^'^tl^ = h'^-'^/''M{L) . (1.21) 



COMPUTER STOCHASTICS IN SCALAR QUANTUM FIELD THEORY 



13 



Fig. 3. A 3D cross section tlirougli a large cluster of a 4D Ising system near the critical point 

We therefore expect finite size scaling of M{L) ~ L'^~^/'^, an anomalous dimension. 
This scaling behaviour should be followed by the clusters, as they are relevant for the 
critical behaviour. They should therefore exhibit a fractal dimension df — d — P/v = 
dj (1 + 1/(5). In 4D we expect mean field critical exponents ^ = \ and thus df — 3. 

In a computer study on a finite lattice the clusters are of limited extension. Small 
clusters are essentially one-dimensional lattice animals, whereas large clusters will 
feel the boundaries and fill the space and then behave like d-dimensional objects. 
Scaling behaviour we can expect only at intermediate size clusters. However, the 
larger the lattice, the larger the domain of intermediate size. In a study of the 4D 
<I>^ theory with 0(4) symmetry [80] the relation between mass s (number of sites 
belonging to the cluster) and size R (the square root of the mean square radius 
between points of the cluster) was studied. 

In fig. H results of simulations for 4D lattices of various size demonstrate the 
existence of a scaling region. Its slope is the fractal exponent in the relation 

s cx R'^i (1.22) 

with a value c?/ ~ 3.05 indeed compatible with the expectations. An earlier study 
for the 3D Ising model [138] had produced a value d/ = 2.51 also compatible with 
the expected relation to the corresponding critical exponents in d = 3. 



2. Triviality of 
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Fig. 4. In this log-log plot [80] of observed combinations (s, R) of clusters at the critical point the 
intermediate scaling region becomes obviously larger for larger lattices. Its slope defines the fractal 
dimension. 

2.1. Rigorous Results 

The claim "This theory is trivial" should be formulated more explicitly: Defining the 
continuum theory as the limit a ^ of the lattice theory, then all n-point functions 
may be expressed through products of 2-point functions (a theory of generalized free 
fields), which even may be those of a free particle (a theory of free fields). 

Let me summarize here the main results. For more details one should refer to a 
review like [51, 20]. 

D < 4: The continuum limit defines an interacting field theory. The renormalized 
4-point coupling is bounded (0 < < const) and perturbation theory 
applies (cf. [49, 50, 15] and further references in [20]). 

D > A: The continuum limit defines a trivial theory with mean field behaviour; 
all connected 2n-point functions (Wightman-distributions) vanish. Most likely 
this is a theory of free fields (cf. [1, 40, 16]). 

£) = 4: No rigorous proof of triviality exists, however there is a collection of results 
which leave little room for another result. 

— Triviality has been conjectured by Wilson [141, 143] in the framework 
of perturbative renormalization group and, under this assumption, many 
scaling results have been derived (cf. [14]). 

— The renormalized 4-point coupling A^ is bounded (0 < A^ < const); for the 
single component theory the critical exponents are bounded from below by 
gaussian (mean field) values [51]. 
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— Considering the equivalence to random walks [132, 133] led to exact in- 
equalities like Xr < 3pis (where pis denotes the intersection probability of 
4D random walks) [1, 2, 3, 40, 28]. 

— Rigorous control of block spin transformations at the critical point (the 
massless theory) proves triviality for small enough Xr [48]; close to the 
critical point there are logarithmic corrections [59, 60]. 

If the continuum limit is trivial and therefore defined by a fixed point of gaussian 
nature then its nature does not change in the complete SU{2) gauge Higgs 
model [62] - although another additional FP at some different value cannot be 
excluded. Assuming triviality one can find regions, where the high temperature 
series (expansion in 1/k) still gives reliable results and the (renormalization 
group) scaling behaviour is already valid. This overlap allows to continue n- 
point functions from the symmetric (hot) phase through the phase transition 
[102]. Triviality implies a bound on the renormalized coupling and on the Higgs 
mass, even if the model is just an effective theory [19, 26]. We will discuss these 
issues below in more detail. 



Fig. 5. Phase diagram of the 1- or (in parentheses) 4-component model 

The phase diagram of the <i>^ model in the (k. A) plane is sketched in fig.^. For 
the 1-component model one has a massive Higgs boson in both phases. At the phase 
transition the dimensionless mass vanishes due to a 0. The 4-component model 
has four degenerate massive bosons in the symmetric phase. In the broken phase 
one finds one massive boson and three massless Goldstone modes (although a > 0). 
The figure also gives regions of applicability of the standard perturbative and non- 
perturbative methods. Monte Carlo calculations cover the whole domain of interest 
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but they are limited to finite volume. On the other hand, the analytic methods 
are limited to certain domains and have to assume dominance of the gaussian fixed 
point (GFP). A thorough analytic understanding, however, helps to analyze the 
Monte Carlo results in a reliable way, as I try to demonstrate below. 

Before we turn to Monte Carlo methods and results let us briefly review the results 
of renormalization group. Perturbation theory gives the /3-function that describes 
the rescaling properties of the coupling with respect to a change of the scale /x. For 
the 1-component model the one-loop result is 

Solution of the differential equation gives 

'^^^^l-A(A)-L.lnWA) '^^^=l + A(M)TB^lnOVA)- ^'''^ 

The cutoff A corresponds to the inverse lattice spacing 1/a, thus A(A) — \bare- For 
the physical scale ^ one could choose some low-energy mass (like the Higgs mass), 
thus A(/i) — Xr- Taking these results at face value, for fixed A(/i) the cutoff A cannot 
be pushed to arbitrary high values due to a zero in the denominator (Landau pole). 
Or, to rephrase this problem as a question: How large do we have to make the 
microscopical coupling Xbare in order to have a non-zero macroscopical coupling Ar? 



Fig. 6. The renormalization of the quartic coupUng as expected from renormalization group 
arguments. 

In fig.^ this situation is summarized. Trajectories of fixed bare coupling are 
bounded from above by the Xbare — oo values of A^. Holding A^ fixed but attempting 
to raise the cutoff (approach the critical point, i.e. the continuum limit) is not 
possible - one gets stuck at the upper bound at some value of A (or, equivalently, 
at some non-zero value of «(/«) > 0. The bound forces the renormalized coupling to 
approach zero in the continuum limit. Only at finite cutoff non- vanishing values of 
Xr are permitted. Thus even a trivial theory may have an effective interaction below 
the cutoff. 
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This discussion is based on the one-loop perturbative result and the renormaliza- 
tion group equations. It was one of the challenges of the Monte Carlo approach to 
confirm (or not confirm) this scenario. Although the asymptotic shape of the upper 
bound is given, the multiplicative factor is of non-perturbative nature. Also, the 
shape at small values of A is regularization dependent. The corresponding ampli- 
tudes are universal up to regularization dependent corrections 0{p? / h?). 

In a sensible regularization the cutoff should be larger than the masses of the 
physical states, thus A//i > 1. If the situation is like discussed above, this provides 
an upper bound on A^. Perturbation theory in the electro-weak gauge coupling g 
relates the renormalized coupling with the ratio between the Higgs mass niH and 
the mass of the heavy vector boson W , 

K = \r\ i?=^ = f^, (2.3) 
2 b 2 mw 

(As mentioned, sometimes the continuum notation with an extra factor 6 is used: 
Ar = 3i?^, cf. table |.) Given the experimental values {F ~ 246 GeV, or mw 
= 80.6(4) GeV, ~ 0.4) the non-perturbative upper bound on the renormalized 
4-point coupling provides an upper bound on the Higgs mass [19, 26], unrelated to 
the usual unitarity-bound. 

What evidence can we expect from Monte Carlo simulations? There are three 
groups of results 

(a) Phase structure and critical exponents, derived from bulk quantities like the 
internal energy, susceptibility or specific heat; detailed investigation of the finite 
size properties of the partition function allows quantitative statements on finite 
size scaling, the thermodynamic limit and logarithmic corrections. 

(b) Monte Carlo renormalization group (MCRG) methods allow to determine the 
flow of coupling and the fixed point structure, including the critical exponents. 

(c) Propagator measurements lead to the mass spectrum and values of the renormal- 
ized couplings. These results are of particular interest concerning the triviality 
bound on the Higgs mass and its regularization dependence. The energy spec- 
trum also allows the determination of phase shifts. 

In the subsequent discussion of this section I first present more qualitative results 
from real space renormalization group, then a high precision determination of the 
partition function and the logarithmic corrections. In sect. 3 the non-perturbative 
results on the bound are discussed. 

2.2. Real Space Renormalization Group 

In perturbative renormalization group the number of couplings considered is very 
small, although the number of couplings generated in the transformation may be 
infinite. The Monte Carlo renormalization group method, also called real space 
renormalization group (RSRG) method since it works in position space, allows to 
deal with a larger coupling space, although one has to truncate as well. 

The Scaling Hypothesis suggests, that the singularities of the thermodynamic 
quantities are due to the divergence of the correlation length [37, 4, 22]. At the 
critical point the system should be scale invariant and under scale changing trans- 
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formations 

T(A,A') 

Z(A) -Z'(A') (2.4) 

the physical properties (i.e. the properties at long distance in terms of lattice spac- 
ing) and thus the partition function should be invariant. Here the transformation 

T is clement of a semigroup. On a lattice it is realized by the Wilson-Kadanoff 
block-spin transformation [143], where the local degrees of freedom are "integrated 
out". We have 

scale transformation A — > A' 
linear factor 6 |A|/|A'| = 

invariance Z{A) = Z'{A') 

field variables $ e T — > $' e T' 

correlation length n^a — > n'^a' — {n^/b)ba 

For lattice systems T transforms an ensemble of configurations on A into an ensemble 
of configurations on A'. The correlation length changes by the scale factor b. 

Given that the original ensemble was distributed according to a Gibbs measure 
for an action S{k, $) let us for the moment assume that the blocked ensemble also 
follows such a measure, but for a different action S{k' , $') (cf. the discussion further 
down). Let us formulate the BST (block spin transformation) $ — > by introducing 
a normalized transformation probability P with the properties 

(a) normalization: / d^A'{'^')P{^' = 1 

(b) positivity: P($', $) > 0; with (a) this implies that for all <I> there is a such 
that P($', $) ^ 0, i.e. there is a non- vanishing probability of mapping into each 
sector of V . 

(c) symmetry: all symmetries of the action S should be respected. 
We then may describe the transformation by 

Z{K) = J d/ZA($)exp[S(K,$)] = J rfMA($)d/XA'(^')^'(^',^)exp[5(K,$)] 

= J diiA' ($') exp[5(«', $')] = Z'{A') (2.5) 

with 

exp[5(K', $')] = J dfiA^Pi^', $) exp[5(/^, $)] . (2.6) 

The BST is largely arbitrary [17] and may be chosen optimal from the point of view 
of simulational simplicity. One has considered both, deterministic and probabilistic 

forms, 

P($',$) = l[S{^'^,,f^>{^)) (deterministic), 

x' (2.7) 
= J|exp($;,,^p«/^'(*)) (probabilistic). 

x' a 

A deterministic transformation with a linear relation between the spins and the 
blockspins makes the new spin at x' an average over the old spins within some block 
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Bx', 

fx'm^ (2-8) 

whereas a nonhnear form could be chosen to keep certain expectation values (in the 
example below: ($^)) fixed during the transformation, 

For Ising spin systems this choice gives the majority rule. Probabilistic transforma- 
tions have adjustable parameters pa that may be used to optimize certain features. 

Under repeated transformations the (dimcnsionless) correlation length shrinks 
and the equivalent action changes. 

/ // * f trivial FP ,„ 

^< ■■• = I oo critical FP ' (^.10) 

k->k'->k"-> ... K* (fixed point) . (2.11) 
The couplings renormalize according 

k' = k'{k) (2.12) 

which serves as the definition of the BST in action space, and, assuming simple FP 
(fixed point) properties, k* ~ k'(k*). If wc start at a critical point, the correlation 
length is infinite all the way: One stays on the critical hypcrsurface. If the correlation 
length was finite it will approach zero in the iteration of the BST, corresponding to 
the approach to a trivial (hot or cold) FP. 

At the FP we may linearize ( 2.12| ) to obtain 



<{^)=<+np{K-K*)0+o{\\K~K*\\) with T„;3=(^|^) (2.13) 

or, for the n-th step 



T., - T^, ■ (2.14) 

The eigenvalues {A} of the FP values of the linearized BST decide on the number 
of relevant (A > 1), marginal (A = 1), or irrelevant (A < 1) operators (couplings) 
at the critical FP. (There may also be redundant operators. They occur e.g. if the 
set of couplings is overcomplete [109, 58]). Fig. |^ shows an example for the RG 
flow. Starting at the FP with a small perturbation, the flow will lead away into the 
direction of relevant operators. This defines the rcnormalizcd trajectory, along which 
there are no contributions of irrelevant operators, i.e. no corrections to scaling. 

Usually almost all couplings (corresponding operators in the action) are irrele- 
vant. The attraction domain of a FP is high-dimensional and only a few (two or 
three) coupling are relevant. This behaviour is the reason for universality classes: 
The same critical surface is domain of attraction for many different initial actions. 
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Fig. 7. Schematic display of the renormalization flow in coupling space; starting from the vertical 
axis at K2 = 



Although the position of the critical points may be different, the critical behaviour 
is not, as it is always given by the properties of T* . 

The linearized BST T* factorizes into even and odd subspaces, depending on the 
symmetry properties of the corresponding terms in the action under reflection $ ^ 
— The relevant eigenvalues are related to the corresponding critical exponent. 
This may be seen from the scaling behaviour of e.g. the correlation length. Assuming 
that its critical behaviour is given by some critical index u we find 

an) = b^'{K'{K)) ^ 
(k - K*)-" = b {k'{k) - K*)'" = bX-'in - n*)'" (2.15) 



(2.16) 



which is consistent only for 

_ ln6 
^~ ha ■ 

If there are marginal eigenvalues there is a line of FPs in the critical surface. 

The position of the FP on the critical surface may depend on details of the BST 
[130]. In particular it may be moved into directions of redundant operators [38]. In 
practice the amount of redundancy is often not clear. Improvement may be desirable 
for both, the FP action (cf. [58]) and the renormalized trajectory [130, 47, 92]. 
Recent attempts to optimize the BST such that the FP action becomes particularly 
simple and as local as possible seem to be promising for asymptotically free theories 
[69]. 

Before I continue, let me mention possible difficulties and caveats. 

— The transformation may introduce spurious singularities as discussed by Grif- 
fiths and Pearce [56, 57]. 

— Some BSTs at first order transitions may produce non-Gibbsian measures ([136, 
137], cf. however [106]). 

— The FP may describe a nonlocal action (no exponential damping of nonlocal 
terms) or it may be at infinity in the parameter space of the action. 

— Depending on the lattice geometry and dimensions only a limited range of block- 
ing factors is available. 
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Like in perturbative renormalization group the number of couplings generated 
in the transformation may be infinite. In practical calculations one therefore 
has to rely on truncations. 
For the following discussion we write the action as a sum of operators, 



S^Yk^So,. (2.17) 



Expectation values and correlators are related to derivatives of the free energy 

= 7^^-^^^P^^{S.Sp)-{S^){Sp) . (2.18) 

BSTs provide different possibilities to obtain information on the critical properties. 

Critical exponents: The cross-correlations between observables measured on the 
original and the blocked configurations allow to get estimators for the truncated 
T-matrix [105, 128]. 

= - (5i"^>(4™^) • (2.20) 

Here T is obtained from ( 2.19| ) by inversion. Due to the truncation there will 
be inherent errors (see however [121]). This approach has been very successful 
for 2D and 3D spin models where it has provided the maybe best estimates of 
critical exponents [18]. 

Flow in the space of observables: Repeated EST leads one to the renormal- 
ized trajectory. Matching of observables allows to relate couplings of different 
simulations. 



Flow in the space of couplings: This would allow to identify FPs and univer- 
sality class in a very intuitive and direct way. 

Unlike in perturbation theory, in non-pcrturbative RSRG calculations it is no 
straightforward task to measure renormalized couplings. It is simple to construct 
the blocked configurations (in position space), but one has no direct means to know 
about the corresponding action. There is the problem to learn the values of (many, 
in the thermodynamic case infinitely many) couplings of the blocked action, which is 
not known directly. A possibility is by operator matching: Matching the observables 
determined by Monte Carlo simulations on lattices of the size of the blocked lattice 
[122, 75, 63, 12]. However, for this one has to combine statistically independently 
determined expectation values which needs high statistics. Another, very efficient 
method has been suggested by Swendsen [129] and amounts to determine conditional 
expectation values in the background of given configurations. 
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MCRG calculation for the one component model have been performed on 
lattices up to size 16* and the results were compatible with the expected scenario: 
one relevant even and one relevant odd direction in operator space. The result for the 
even critical exponent was good {ly = 0.53(4) from e.g. [11, 90], or from [21]) but the 
statistical accuracy not sufficient to identify the expected next-to-leading marginal 
eigenvalue. No convincing signal for logarithmic corrections could be determined in 
those calculations. 

However, from the flow structure of the renormalized couplings one could try to 
find out, whether different values of the bare quartic coupling are in one universality 
class. This was done in [90] in a space 24 different couplings corresponding to 6 
quadratic (local and distance 1, 2, 3, and 4 spin-spin products) and 18 quartic (local 
and neighbouring 2-, 3- and 4-spin products) operators. The simulation was done 
on a 16* lattice with most of the initial couplings set to zero. Only the nearest- 
neighbour coupling K and the bare 4-point coupling A were non-zero. For given fixed 
values of A (0.2, 0.4 and 2) the model was simulated for the critical value of k and 
BSTs were performed down to lattices of size 8* and 4*. The renormalized couplings 
in the 24D space were determined with Swendsen's method [129]. 



Fig. 8. Flow structure in coupling space; starting on the critical surface in the (ki,K2) plane 
(corresponding to the bare quartic coupling and n.n. coupling) the RG flow moves towards a 
common FP, independent of the bare quartic coupling; it is irrelevant. 

The observed flow may be summarized as in fig. |^. From different starting values 
of the bare quartic coupling on the critical surface the flow is towards a common 
FP. In the study different BSTs were considered. As the FP position may depend 
on the particular BST one could expect different flow structure, which indeed was 
observed [21, 89, 90]. The deterministic nonlinear type of BST discussed above 
produced a flow towards the FP characteristic for the Ising system. The overall 
result demonstrated, that the bare quartic coupling is not relevant for the FP and 
thus for the continuum theory - in agreement with the notion of triviality. 

2.2.1. Critical Indices from Bulk Quantities 

Bulk quantities like the internal energy per link {E) or the magnetization ($) are 
the most direct way to the determination of scaling behaviour. However, they are 
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affected by finite size effects and therefore the critical exponents may be masked 
behind the smooth crossover between the phases typical for finite systems. Finite 
size scaling [10, 4] relates finite size properties with the critical exponents in the 
thermodynamic limit. However, usually one needs very precise determination of 
the quantities involved. The histogram method, an old idea [30] reshaped into 
the multi-histogram approach by Ferrenberg and Swendsen [31, 32, 33], allows to 
combine measurements at different couplings in an efhcient way. 

We rewrite the partition function in terms of the density of states (for simplicity 
we assume that the set of energy values is discrete, otherwise we have to rely on 
some binning of the allowed range). 

Z{k) = J d/XA($)exp[KS'($)] = J dn/^{^)J2S{S~E)exp[KS{^)] 

E 

= P{E) eM^E) = ^p(i?|/t) . (2.21) 

E E 

The probability for a configuration having the total action E for coupling n is p{E\n) 
and is estimated by the corresponding histogram values /i(k, E) obtained in the 
Monte Carlo simulation of the system. Simulations at different values of the coupling 
provide different estimators for the spectral density p{E): 

Po{E) = h{Ko,E)exp{-KQE) , 
Hi pi{E) ^ h{Ki,E)exp{-KiE) , etc. (2.22) 

However, the statistical accuracy for these different estimators will be best for those 
values of E close to the corresponding expectation value (i.e. near the peak of 
the corresponding histogram). A suitable combination of histograms [31, 32, 33], 
however, gives a reliable estimator over a broader range of energies (and therefore 
couplings). From this estimator for the spectral density one may construct energy 
histograms at arbitrary values of the coupling, as in fig. |^. 

This technique may be generalized to other observables and terms in the action. 
Taking into account energy and magnetization we have the partition function 

Z{k, h) =Y^ p{E, M) exp{KE + hM) (2.23) 

E,M 

allowing the precise determination of expectation values of arbitrary functions of E 
and M, 

{f{E, M)) = \Y. P^^^ ^^)/(^' exp('^i^ + hM) . (2.24) 

E,M 

An example is the specific heat Cy ; in fig- ^ both, the values determined at specific 
couplings and the continuous curve from the multi-histogram analysis is shown. Each 
curve contains all the information of the corresponding individual couplings (for that 
lattice size). 



The representation (2.23) also allows for analytic continuation to complex val- 
ues of K and h since only the exponential has to be continued. One therefore can 
determine the zeroes of Z in the complex k and h planes. In the actual calculation 
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Fig. 9. The logarithm of the probabiUty distribution (histogram) as reconstructed from the 
multihistogram: ln[p(£;)e''o^^] vs. E, for kq = 0.15055 and lattice size 16^. 



Fig. 10. The specific heat Cv for the one-component model determined for lattice sizes 8* 
12*, le'', 20*, and 24* (from [84]). 



this can be done reliably only close to the values of the couplings, where the original 
histograms have been determined, otherwise the statistical errors of the estimated 
spectral density blow up and lead to large uncertainties. Experience shows, that ex- 
trapolation is feasible in an elliptical domain with the real positions of measurements 
between the foci of the principal axis. 

The study of zeroes of the partition function has been started by Lee and Yang 
[146, 147]; zeroes in the complex external field h are called Lee- Yang zeroes. They 
have to lie at purely imaginary values of h (on the unit circle in the fugacity z = e~^^ 
plane) [146, 147]. Zeroes in the even coupling (the temperature) have been studied 
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in particular by Fisher [36] and are thus usually called Fisher zeroes. On a finite 
lattice the partition function is a finite-degree polynomial in e"^** and z. 



Fig. 11. Left: Finite size scaling behaviour of the imaginary part of the closest Fisher-zero k\ 
for lattice size L = 8, 12, 16, 20 and 24 at /i = 0. (a) The leading exponent is related to the 
slope of the straight line and a fit gives -2.088(6) slightly differing from the mean field expectation 
— Xjv = —2. This difference is due to the presence of logarithmic corrections as demonstrated in 
(b) where the leading behaviour has been cancelled by considering L^Im/ti . Right: Finite size 
scaling behaviour of the imaginary part of the closest Lee- Yang-zero hi for lattice size L = 8, 12, 
16, 20 and 24 and k = 0.149703. (a) The leading exponent is related to the slope of the straight 
line and a fit gives -3.083(4) slightly differing from the mean field expectation —1/5 = —3. This 
difference is due to the presence of logarithmic corrections as demonstrated in (b) where the leading 
behaviour has been cancelled by considering L'^lmhi (figures from [86]). 

Expressing, e.g. Z through a product of zeroes we find, that the specific heat - 
the second derivative of its logarithm - is peaked for those real values of k that are 
closest to a nearby complex zero: 



In the thermodynamic limit the zeroes accumulate and pinch the real K-axis (and, 
for critical coupling, the real point h = 0): They build the critical singularity. 
The A-dependence of the positions of the zeroes closest to the real coupling axis 
can be related to the thermodynamic critical exponents (finite size scaling) and the 
logarithmic corrections. 




(2.25) 
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Assuming that the gaussian FP controls the leading critical behaviour, perturba- 
tive renormalization group methods may be applied to predict the critical properties 
[14] and the finite size scaling behaviour of the partition function zeroes [13]. For the 
one-component model this was done recently [84, 86, 85] (for the 0{N) model 
cf. [83]) with the following results for the Fisher zeroes (using the reduced coupling 
t) and the Lee- Yang zeroes: 

\tj\ - L-2(lnL)-3 , (2.26) 
\hj\ - L-^{\iiL)-i . (2.27) 

Here the leading exponents are the negative inverse of the mean field values for i/ and 
(5, the multiplicative logarithmic corrections depend on the number of components of 
$ and have other values for the A^-component model. This result also gives the finite 
size scaling behaviour for the peak values of the specific heat and the susceptibility 

CLit = h = 0) - (lnL)3 , (2.28) 
XLit ^h^O) - L^ilnL)^ . (2.29) 

In a high statistics simulation the multi-histogram techniques were used to study 
bulk quantities of the model and the positions of the closest partition function 
zeroes [84, 86, 85]. The numerical result together with the Ferrenberg-Swendsen 
technique allowed sufficient precision to identify both, the leading critical exponent 
and the form of the logarithmic correction. The results for the Lee- Yang and for 
the Fisher-zeroes (cf. fig. 0) confirmed the expectation. The leading even and odd 
exponents agree with the mean field r esult. The logarithmic corrections are in good 
agreement with the expressions ( ^.26 ) and ( ^.27| ) derived on ground of perturbative 



renormalization group. 
3. Consequences of Triviality 

Let us now concentrate on results derivable under the impact of triviality like the 
bound on the renormalized coupling and the Higgs mass. We therefore discuss in 
this section in particular the 0(4) model because of its close relationship to the 
unified theory of electro-weak interaction. The notation refers to the lattice action 
(O) given above. 



3.1. Results from Non-Monte Carlo Techniques 

Using a local potential approximation and Wilson recursion relation technique for 
renormalization group calculations Hascnfratz and Nagcr [68] obtained analytic and 
non-perturbative results on the bound for the one-component model. The effect of 
the inherent approximations is hard to judge in that approach. A real breakthrough, 
however, were the results of Liischer and Weisz [100, 101, 102] and for this reason I 
want to briefly summarize them. 

The authors study the A^(> l)-component model on a hypercubic lattice. Start- 
ing from the symmetric phase they proceed as follows (cf. fig. 

Symmetric phase: The high temperature (small k) expansion for the renormal- 
ized mass m^, coupling A,, and the critical value Kc is determined up to 0{k^^). The 
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values seem to be reliable up to relatively small values of ~ 0.5, already within 
the scaling region. 

Scaling (critical) region: Here one uses the results of perturbation theory in the 
renormalized quantities: rrir and for fixed bare A, using 3-loop /3-functions. 

One has to assume that the gaussian FP defines the universality class. The a priori 
unknown integration constants like Ci in 

n<K,: am,(A) = (A)(/3iA,)-^^/'3i exp (^-^) {1 + 0(A,)} (3.1) 

contain the non-perturbative information and are taken from the K-expansion values 
in the overlap region. Corrections to scaling are 0(a^ (In a)'^). The scaling form is 
used to continued through the phase transition to the phase with broken symmetry 
in a region < 0.5. 

Broken symmetry phase: The Integration constants below and above the critical 
point are related through the massless theory at the critical point, 

K> Kc:C[{X)=e^Ci{\) . (3.2) 

They are renormalization prescription dependent. 



Fig. 12. Quantitative plot of renormalization group trajectories for the 0(4) model (from [102]). 
The abscissa variable A corresponds to values of the bare coupling < A < oo . The range of k- values 
roughly covers the scaling region \mr\ < 0.5. The curves correspond to constant renormalized 
4-point coupling — they approach the critical line, but they never touch it; In the continuum limit 
one is forced to vanishing values of A,.. 

This program exploits standard analytic techniques in a very efficient way and 
worked nicely [100, 101, 102]. In fig. |l^ trajectories of constant renormalized 4- 
point coupling are shown: They never cross the critical line, therefore only a non- 
interacting continuum limit is feasible. 

We summarize the restrictions or assumptions for this approach: 
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— the K-series is truncated; 

— triviahty (gaussian FP) is assumed; 

— vahdity of perturbation theory in the renormahzed couphngs (renormahzation 
group at the 3-loop level) is assumed. 

Direct Monte Carlo simulations can check these assumptions and confirm the validity 
of the approximations. They also can determine further results concerning, for 
instance, the regularization dependence. A consistent picture gives weight to both, 
the analytical and the numerical approach. 



3.2. Monte Carlo Simulations 

First results go back to 1984 [140, 135] but stayed widely unnoticed. Only with the 
very high statistics studies of various groups [65, 108, 87, 88, 7, 74] since 1987 the 
results became convincing enough. These studies worked with lattice sizes up to 
16^. Most results have been obtained in the broken phase, both for j = and j > 0. 
Typically several 100.000 sweeps and measurements per point have been performed, 
using first the Metropolis, later the cluster algorithm. In summary thousands of 
hours supercomputer CPU-time (like e.g. Cray-YMP or equivalent) have been spent 
for these analyses. Some studies were concerned with the one-component model but, 
again, most concentrated on the 0(4) model. 

For the determination of the rcnormalized coupling in the phase of broken sym- 



metry we need (cf. (2.3) and table ||) 

— The vacuum expectation value E = (1^9) ; the complication is, that for j — 
in finite volumes there is no spontaneous symmetry breaking and one has to 
consider suitable definitions. 

— The field renormahzation constant Z or, equivalently, the pion decay constant 
F = E/a/Z. This may be determined e.g. from the propagators of the massless 
Goldstone modes. 

— The rcnormalized mass of the massive (Higgs) mode, determined from the 
measured Higgs propagator. 

In addition to these quantities necessary for the derivation of A^. we may also want 
to identify clear signals for the existence of Goldstone modes, like e.g. propagators 
at J = and j > 0, just to be sure that the symmetry breaking mechanism is the 
conjectured one. 

The problem to define, what one means by spontaneous symmetry on finite lat- 
tices, turned out to be hard in practice. A straightforward definition originally 
employed was 

1 ^ 



S = (|(^«|) where ^« = — ^y,^, (3.3) 



and rUr, F and Z were determined from the propagators of the field components 
parallel and perpendicular to the value of the corresponding configuration, which 
meandered around in the 0{N — 1) subspace. That seemed to work nicely. In fig. 

we show results for the Higgs mass determined that way on lattices of various 
size. Further results showed that the field renormahzation constant stays close to 
the value 0.97 from the PT up to raj. k, 1. 
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Fig. 13. The dimensionless renormalized Higgs mass a(K)mr determined on lattices of various 
size and extrapolated to infinite lattice size (from [113]). 



A better understanding of the infinite volume limit of the ad hoc ansatz (3.3) 
was important to achieve systematically improved quality of the results for S. Now, 
for J ^ the symmetry is broken explicitly, all quantities are well defined and can 
be measured without reference to spontaneous symmetry breaking. For small values 
of the external field the would-be Goldstone states are still very light and therefore 
they feel the finiteness of the volume best. We have a situation of two different 
kinds of long correlation lengths (fig. |l^). At the 2nd order PT a{K)mr — > 0, 
because a(K) — > for k — + Kc, since the theory is critical. Along the 1st order PT 
a{K)mG{j) because the Goldstone boson mass mcij) for j 0. In the 
discussion below we will denote masses of the massive mode mg. and of the Goldstone 
mode TOtt in analogy to the nonlinear cr-model. 



Fig. 14. Phase structure in the (k, j)-plane. 



Let US recall again the important parameters of the model. The field has four 
real components. In the symmetric (high temperature) phase the particle spectrum 
has 4 degenerate states. In the low temperature phase k > Kc, at infinite volume 
0(4) symmetry is spontaneously broken. In this phase 



(^2> = S, (^^)=0 (* = 1,2,3) 



(3.4) 
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Here the spectrum of the theory contains three massless Goldstone bosons corre- 
sponding to the excitations in the 0(4) directions a = i = 1,2,3. Therefore, the 
correlation functions do not fall off exponentially at large distances, but with an 
inverse power of the distance. Specifically, the Goldstone boson two-point function 
satisfies 

lim in^\x-y\^{v>ivi)=Z6^'- (3-5) 

\x-y\—^oo 

For 0(4) symmetry the model contains six currents, all conserved at j = 0, whose 
charges generate the group. The corresponding Ward identities strongly constrain 
the behaviour of the correlation functions at large distances. In fact, the asymptotic 
behaviour of all the Green functions associated with the currents and with the fields 
(fi'^ is determined by the two low energy constants E and F [41]. In analogy with 
QCD, F is often referred to as the pion decay constant. It specifies the matrix 
elements of the axial currents Ajj^{x) (in Minkowski space) between the ground state 
and single Goldstone boson states, 

(0|4(0)|7r'=(p)) = i6^%F . (3.6) 

The low energy constants Z, S and F are related through 

F = Y./v^. (3.7) 

For j > the Goldstone boson mass m,r has a non- vanishing value which in the 
lowest order of chiral perturbation theory is 

ml =3^. (3.8) 

Whenever (almost) massless states dominate the system the low energy (long 
distance) structure is of a universal character determined by the symmetries of the 
dynamics in terms of a few low energy constants and by the geometry of the finite 
system. One therefore may study an eff'ective field theory with only Goldstone 
modes and the symmetry structure of the underlying model 0{N) 0{N — 1) 
[95, 43, 67, 111]. The effective Lagrangian has the same symmetry as the original 
system but is simple enough to permit a systematic perturbative expansion. It has 
the form 

^eff = \F^{d,Sd^S)-^HS 

+further terms of higher order in H and S , (3-9) 
Z{H) = j dij{S)exp{ j dxL^s) . (3.10) 

In higher orders of this perturbative Lagrangian some further low energy constants 
Ae, Ai? and Am are required [67]. They are the scale parameters determining the 
logarithmic dependence of S, F and on j. In systems of finite size L this 
expansion amounts to an expansion in powers of L~^. It is reliable because the 
interaction between the Goldstone modes is weak at low momenta. 
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The first application of this idea [42, 95] has been made in the context of the 
study of finite size effects in QCD with light quarks, whose low energy properties are 
determined by the chiral SU(2)(g)SU(2)~0(4) symmetry [41]. Therefore one refers 
to this approach as " chiral perturbation theory" . A comprehensive description can 
be found in [43, 94, 67]. 

The effective theory and its expansion relates values of quantities determined at 
finite volumes and j > with those at j = at infinite volume, 

Chiral Pert. Th. 

E(j, V), F{j, V),Z{j, V) S, F, Z. (3.11) 



Fig. 15. Tile typical domains of application for two type of expansions indicated by lines of 
constant m^lrriT, in the (k, j) plane. The dashed curve corresponds to m^L = 1 on a 10* lattice. 
The regions below (above) the dashed curve are the U (V) domains for L = 10 (from [64]). 

There are two types of expansions, valid in different domains with a non- vanishing 
overlap (cf. fig. 

{JX) TOtt ^ 1/i : This domain is characterized by very small symmetry breaking, 
where the symmetry is restored. The finite size effects are large here. The 
expansion is in powers of keeping the total magnetic energy 

^^0 = SjL^ , (3.12) 

fixed. Thus j is treated as a quantity of the order 0{L^'^). In 14 the correlation 
length grows as for L — > cxd. 

(V) m,Tr ^ 1/L : In this region j is still small, although larger than in U. Finite 
size effects are smaller than inU. A typical range observed is ma 
Here the symmetry is not restored. The expansion in powers of L^^ in that 
domain keeps fixed 

vo = ^jLyF^ =mlL\ (3.13) 
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For large uq the results of the expansion U should smoothly go over into those of 
the expansion V. There is an overlap of the regions of validity of both expansions 
around ttLt^L ~ 1 [64]. In fig. |l^ the dashed curve indicates where the U and V 
regions in the n-j plane meet for L = 10. 

We do not want to go into details of the various expansions. Just in order to catch 
the spirit of the approach let us give an example for one result of the ^//-expansion. 

where the second term is 0(L^^) and 

The second term here is 0{L^^). The constant f3i = 0.14046 depends on the lattice 
geometry, as does the constant in p = 0.01523 + 0.0095 ©(InLAs). The terms 
0(L~^) are small and involve higher orders in the chiral expansion. 



Fig. 16. Example iocthe type of data obtained for ((/J^ ) (from [64]). Tlie solid curves correspond 
to a fit according (3.14) up to 0(L~^), the resulting value of tiie only parameter E is indicated by 
the full dot at j = 0. 



Eq. ( |3.14 ) gives the functional form of the j and L dependence in terms of the 



thermodynamic quantities S and F. Fig. |1^ gives an example of the type of data 
and fit quality obtained. The leading term in ( |3.14 ) provides the value of S. The 



non-leading term in principle can give F, but is statistically not sufficiently reliable. 
For this parameter one better uses the t-dependence of the Goldstone propagator. 

The values obtained in this way prove to be very precise, as demonstrated 
in fig. |l^. One should keep in mind, that these are the thermodynamic valuesl The 
fit to the scaling behaviour gives [64] 



= 0.672(6)(-i)|ln(-t)|i and = 0.3036(10). 



(3.16) 
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Fig. 17. The values of S obtained from the fit of j > 0, finite volume data to the chiral perturbation 
expression; the curve shows the fit to the expected scaUng behaviour (from [64]). 



TABLE II 

Results on various quantities of the 0(4) <E>^ model as obtained with the action < ]1.3j) in various 
studies. 



K 


S 




F 


Z 


■mo- 


nic/F 


Ref. 


0.3130(7) 


0.196(1) 




0.199(5) 


0.976(6) 


0.50 


2.517(65) 


[102] 


0.3101(5) 


0.166(6) 




0.168(5) 


0.974(7) 


0.40 


2.380(69) 


[102] 


0.3077(3) 


0.130(3) 




0.132(3) 


0.973(5) 


0.30 


2.280(58) 


[102] 


0.3058(2) 


0.093(2) 




0.094(2) 


0.972(5) 


0.20 


2.129(46) 


[102] 


0.3046(1) 


0.051(1) 




0.052(1) 


0.971(5) 


0.10 


1.941(34) 


[102] 


0.355 


0.4036(9) 


0.402 


0.4109(13) 


0.965(2) 


1.09(10) 


2.65(25) 


[64] 


0.330 


0.3027(11) 


0.301 


0.3075(16) 


0.969(3) 


0.92(7) 


2.99(24) 


[64] 


0.325 


0.2769(12) 


0.276 


0.2807(19) 


0.973(5) 


0.81(2) 


2.89(9) 


[64] 


0.310 


0.1643(63) 


0.163 


0.1668(73) 


0.97(1) 


0.39(2) 


2.34(22) 


[64] 


0.3075 


0.132(13) 


0.128 


0.135(14) 


0.96(1) 


0.29(1) 


2.15(30) 


[64] 


0.3100 


0.15549(10) 




0.1584(7) 


0.964(7) 


0.4073(7) 


2.572(15) 


[54] 


0.3080 


0.12827(6) 




0.1301(3) 


0.972(3) 


0.3223(12) 


2.477(14) 


[54] 


0.3060 


0.09045(10) 




0.0915(2) 


0.977(3) 


0.2072(7) 


2.264(14) 


[54] 



In table ^ the values of E obtained by the discussed chiral perturbation expansion 
method [64] is compared with the original "brute force" ansatz (O) and one finds 
surprising agreement - although the new method definitely is of higher quality. 

Let us summarize the different ways to determine the low energy parameters. 

Determination of E: (a) from the j > data with help of the discussed 
chiral perturbation expansion in j and L^^; (b) at j=0 from the distribution of the 
"local" mean magnetization (|^"|) defined in (|3.3D - for a detailed discussion of this 
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constraint effective potential approach cf. [52, 55, 53, 54]. 

Determination of Z: from tlie Goldstone propagator (a) at > from chiral 
perturbation expansion formulae [64] and (b) at j = using the field components 
transversal to the direction defined by the magnetization (|^"|). 

Determination of nia : (a) from the component of the field parallel to the 
direction of (for each configuration) ; (b) from the susceptibility 

X=|A|((|n')-(l^"l)') , (3.17) 

which is in leading order proportional to The latter method appears to work 

very efficiently [54] . 



Fig. 18. Comparison of results for the renormalized coupling for different values of the cutoff 
A = Xliarria) (from [113] in 1990). 

In table ^ we compare the results derived over the years from various groups. 
One observes increasing quality of the numbers. Comparing the Monte Carlo results 
for the renormalized coupling at values mo- < 0.5 with those of [100, 101, 102] one 
finds good agreement (within the statistical errors). This was also confirmed by 
further Monte Carlo results of other groups [87, 88, 7]. The resulting estimates for 



the upper bound on the renormalized coupling are shown in fig. 18 and fig. |19|. One 
clearly notices the progress in the statistical accuracy. 

If we assume that the cutoff (given in units of in the fig ure) should not be 
below ^ 2 then we read off an upper bound oiRiL 2.6(3). With (2_^) this translates 



into VTLH ^i670(80) GeV. This bound comes from the field theoretic properties of the 
$^ model as a direct consequence of triviality and should not be mixed up with the 
so-called unitarity bound on the Higgs mass. The bound expresses the fact, that 
for this part of the standard model we cannot perform the continuum limit. Or, 
put differently, for a given Higgs mass there is a limited domain of validity of the 
regularized theory. If the Higgs mass comes out to be much smaller, e.g. of 
the order of 200 GeV, the bound allows a cutoff of the order of the Planck length. 
Up to the allowed cutoff we still can consider the theory a model with an effective 
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Fig. 19. Results for the renormalized coupling for different values of the cutoff A = l/{am(j). The 
4 points with larger error bars are from [102] and correspond to twice the error quoted, following 
their discussion of uncertainties. The upper group of points corresponds to a different action as 
will be discussed below (from [54] in 1992). 

interaction. However, we conclude that the theory cannot be a truly fundamental one 
(we cannot remove the cutoflF completely) and we therefore may expect constituent 
structure. 

All these results have been obtained away from the critical point, therefore one 
does have contributions from non-leading operators. These corrections to scaling are 
regularization dependent and O(a^). The cutoff then represents new, possibly more 
fundamental physics which would replace the standard-model at higher energies. 
However, its influence is effective on all scales and manifests itself in the cutoff 
dependence of all obscrvablcs. In describing the experiments through a lattice model, 
the higher in energy we go the stronger we expect such cutoff effects to be. More 
and more terms would have to be included in the effective action in order to explain 
experiments with sufficient accuracy. Within the scaling region one expects to find 
these dependencies to be only weak, but this knowledge is mainly based on rough 
estimates and there are only few explicit studies [7, 8, 70, 53, 54]. 

Value and shape of the bound depend on the adopted regularization scheme, i.e. 

on: 

— The lattice geometry; an investigation of the scalar field theory on various types 
of lattices is required [7, 110, 112, 70]. Methods for comparing the results 
obtained by means of different regularization schemes are being developed [66, 
102, 91]; they require calculations of high precision in a region of the coupling 
parameter space where the universal asymptotic scaling associated with the 
("trivial") gaussian fixed point sets in. 

— Additional terms in the action; an example is the introduction of a straight 
distance 2 coupling in order to remove the O(p^) dependence in the lattice 
perturbative propagator (Symanzik improvement). 
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— Other terms of 0(1/ A^) in the action; Heller et al. [72, 73, 74] have studied these 
effects in the 0(A^)-niodel by systematically adding operators of increasing order 
in 1/A, both, in continuum and on the lattice, in particular on an _F4-lattice 
[110]. 

— Other regularization schemes: Recent work demonstrates the regularization 
dependence in the continuum 0(4) model with a higher derivative regulator 
[79, 78]. 

Comparing results from different actions obtained at finite cutoff is not straight- 
forward, though. Different actions have different scaling violation properties. Choos- 
ing couplings such that the observed correlations lengths agree is not the only and 
not necessarily the best definition. Consider e.g. an effective action for block vari- 
ables (which usually is more complicated). Comparing the results of the original and 
the block action at a point in their respective coupling constant spaces where the 
dimensionless correlation lengths agree, the block action will definitely have better 
continuum properties, it has less corrections to the leading scaling behaviour, ft is 
"closer to the continuum limit" . 

In an asymptotically free theory that issue can be resolved by comparing the 
behaviour of the /3-function in the weak coupling region [61, 100, 66]. In our case 
in principle we can rely on renormalized perturbation theory. In fact the propor- 
tionality constant C in the two-loop expansion for the renormalized mass in ( |3.l[ ) 
is regularization scheme dependent. We therefore add a subscript s denoting the 
regularization scheme (lattice action or geometry) . Asymptotically, at small enough 
small Ar, as/o-s' = Cs/Cs' can be determined by lattice perturbation theory. 

For a comparison of different effective theories at finite cutoff and not too small A^. 
we can no longer rely on that approach. It would be appropriate to choose another 
physical observable O, that may serve as a means to measure the distance from the 
continuum limit independently. Comparing results due to different actions Sg (where 
the index distinguishes between the different actions and regularization schemes) at 
the same values of O then provides an estimate for the regularization dependence. 
This observable could be for example a TTTr-scattering cross-section, the cr-decay 
width or the violation of euclidean rotational invariance of some quantities [102, 8, 
70, 54, 93]. Different operators O will usually have different cutoff dependence and 
their choice needs further motivation: There is no "best" observable. Quantities like 
the ratio R should be considered versus such a "distance" O to the continuum limit. 
Of course we could still use the inverse correlation length TOphygCis, but in order to 
compare different schemes we have to rescale at each point such that O agrees. 

In a Monte Carlo study of the 1-component model [93] the violation of eu- 
clidean rotational invariance (in real space and in momentum space) of the Higgs- 
propagator was used as measure for the "distance to the continuum limit" . Four 
kinds of lattice actions, respectively lattice geometries, were considered. 

C : The standard action on the hypercubic lattice (p^). 

C : An action including diagonal nearest neighbour terms with the same coupling 
[7]. 

: The nearest neighbour action on an i^4-lattice [8, 9]. It's embedding into a 
hypercubic lattice may be imagined by removing all odd sites of the grid. The 
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symmetries for this case forbid terms like ^^^pfi in the inverse momentum 
space propagator and guarantee Lorentz-invariance to a higher order in the 
momentum cutoff 0(A~*) than the first two hypercubic lattice schemes which 
had violating terms 0(A^^). 

S : The Symanzik improved action which includes anti-ferromagnetic couplings to 
the next-to-nearest neighbours along the 4 euclidean axis-directions. These 
additional terms arc incorporated in order to remove the undesired 0(A^^) 
contributions. The remaining corrections of 0(A~*) in the lattice dispersion 
relation guarantee tree-level improvement [134]. 

All four lattices may be embedded into the regular hypercubic one. 



Fig. 20. Fig. from [93], showing the amount of rotational symmetry in the inverse boson propagator 
in momentum space, determined for 4 different actions. 

Fig. compares the rotational invariance properties for the boson propagator 
in momentum space. The deviation from the exact continuum symmetry may be 
quantified by considered the coefficients in an expansion in spherical harmonics. Fig. 
^ exhibits the final results for R = rria/F (proportional to the square root of the 
renormalizcd coupling) as functions of a measure for the deviation from rotational 
symmetry. As discussed above, this may be a more appropriate measure than the 
correlation length itself. One find a sizable difference, about 30-40%, between the 
values for the actions C and C (which violate rotational invariance to a larger 
extent) on one hand and and S on the other (which have no 0{p'^) violation in 
the tree- level propagator). 
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Fig. 21. The ratio (2.3) vs. a measure of rotational invariance for the one-component model; the 
curves are fits to the asymptotic behaviour; the left-hand part give the results on an extended scale. 
The symbols used denote the different action (C: bars, C": squares, F4: triangles, S: asterisks) 
(fig. from [93]). 



In the results for the 0(4) model for [8] and 5' [54] the authors estimated the 
regularization dependence by comparing the results at given value of another mea- 
sure of rotational symmetry, derived from a perturbation expansion of the Goldstone- 
Goldstone scattering cross-section. There the difference between F^, S compared to 
C was smaller than that observed in [93], of the order of 10-20 % for 1-10% violation 
of rotational symmetry. Although this variation could be due to the study of differ- 
ent models (one component vs. four component fields) it may also indicate a strong 
dependence on the choice of measure for rotational invariance. Maybe systematic 
studies of higher order terms, like those initiated recently [74, 72, 73, 79, 78] are 
helpful to clarify this issue. 

Let us conclude this chapter with a summary of the Monte Carlo results for the 
0(4) $4 model. 

— The observed scaling behaviour is consistent with the triviality conjecture. 

— One clearly identifies Goldstone bosons and the massive Higgs mode. 

— One finds good agreement with analytically derived results in regions, where 
the comparison is possible {aniH < 0.5). 

— Finite size scaling works. 

— Chiral perturbation theory provides an excellent description of j > finite vol- 
ume results and allows a reliable extrapolation to the thermodynamic situation 
of spontaneous symmetry breaking. 

— Allowing for a cutoff > 2mH (or "continuum-like behaviour" up to a few per- 
cent) one obtains an upper bound for the renormalized coupling, that translates 
(using the experimental values for the gauge coupling and the VF-mass) into the 
bound on the Higgs mass [71], 



TJiH < 710(60)GeV. 



(3.18) 
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— Some (relatively weak) regularization dependence of this bound (and other 
quantities) has been observed. 

— Further work includes fermions (Yukawa models), but no drastic change of the 
bound has been found. 

The recent review of Heller [71] summarizes the situation and gives further rele- 
vant references. 

4. Scattering: Phase Shifts and Resonances 

Most observables that have been considered in lattice calculations so far are masses 
(from propagators) and coupling constants. Most experimental measurements in 
high energy physics, however, lead to cross sections and their analyses to scattering 
phase shifts. Indeed, resonances are described by the slope and position of the values 
of the phase shift passing through 7r/2 and the width is a quantity derived under 
certain additional assumptions (like the Breit-Wigncr shape). In fact, considering 
the 0(4) model we are in a situation, where the Higgs boson is a resonance 
(coupling to Goldstone states) and, strictly speaking, no asymptotic state of the 
theory. For the determination of the Higgs mass it is therefore not sufficient to just 
determine the single particle propagator over larger and larger distances. In that 
limit one will only learn about the lowest lying energy level in that channel, which 
is not the Higgs mass but a 2-particle threshold. 

In finite systems the energy spectrum is discrete and may be related to the 
phase shifts in the infinite volume system [96, 97, 103, 99, 98]. The leading volume 
dependence is connected to the scattering lengths and early calculations in the 4D 
Ising model determined the s-wave scattering length [108, 39]. Meanwhile there are 
several studies: 

— In a pioneering work Liischer and Wolff obtained phase shifts in the 2D 0(3) 
nonlinear cr-model [103]. There the result is known analytically from conformal 
field theory and one found agreement with the non-perturbative calculations. 

— Below I discuss another 2D study [45, 46, 44] with the aim to determine phase 
shifts in a resonating situation and to find the resonance parameters in a non- 
perturbative way. 

— Fiebig et al. [34, 35] investigate the 3D QED and determined phase shifts for 
the scattering of mesonic bound states (in the quenched situation, i.e. without 
dynamical fermions in the vacuum). 

— In the 4D 0(4) model there have been studies in the symmetric phase [114] 
and in the broken phase with an additional external symmetry breaking field 
[148, 149, 139]. The latter model is a test case for resonating phase shifts as well, 
since one has a situation with light particles (the Goldstone bosons have non- 
zero mass due to the external field) and one heavy resonance (the Higgs particle 
decays into the Goldstone states and is no asymptotic state of the theory). 

Only the efficient algorithms available for bosonic systems allow the high precision 
that is necessary to determine the energy levels and phase shifts, as discussed below. 
For this reason the scalar field theories are the leaders in that approach - the reason 
to discuss approach and result in this last section of this review. Also, in that process 
a lot can be learned on the representation of physical states an the lattice. 
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(a) 



(b) 



© 



© 



^ ^ ^ 



Fig. 22. The spectral density in the s = VF^-plane; (a) infinite volume, (b) finite volume. 



4.1. Energy Eigenspectrum of Finite Systems 

Consider the propagation of a state \N) in a channel with 1-particle (mass M) and 
2-particle states (mass m = 1, i.e. we choose the hght state as the mass unit for 
simphcity of notation, where not stated otherwise). In infinite volume the analytical 
structure of the propagator in the relativistic variable s is as shown in fig. |2^(a) , with 
an elastic cut starting at 4to^ = 4, higher inelastic thresholds and (for M < 2m) 
a bound state below the elastic threshold. The dispersion relation representation 
gives 

T{s) = po\c+ - r ds' (4.1) 
or, Fourier transformed to euclidean time t, 

Git) = CQg-*" + - / dW'e^^ 'p{W') (4.2) 

71" J 4 

where W denotes the energy and p is the spectral density. If the single particle 
state is above threshold it describes a resonance and is no longer part of the discrete 
spectrum and no longer an asymptotic state. 

In finite volume the continuous spectrum is replaced by a discrete set (cf. fig. 

Git) = coe-^'' + J2 c^e-'^^\ (4.3) 

i 

with the eigenspectrum {M^i}. A possible resonance is represented by the form of 
the discrete spectrum and not by a single energy level like a bound state. 

How can we determine the eigenspectrum in realistic lattice calculations? The 
answer is given by standard methods of quantum mechanics. One has to find a 
(sufficiently) complete set of operators iV„ with the correct quantum numbers (sym- 
metry properties) of that channel. Each iV„(i), n = 1, 2... is defined on a single time 
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slice, and t denotes the separation of the time shoes. In the Monte Carlo simulation 
one then measures the connected correlator, 

Cnm(t) = (iV„*(0)iV™(i))c = {Nn*{0)Nrn{t)) " (iV„*)(iV™) . (4.4) 

The transfer matrix formalism (see e.g. [77]) yields the spectral decomposition 

Cnm{t) = J2{K\l)e~'^'Hl\N^)c (4.5) 
= ^.^^.We-'^'*. (4.6) 

For simplicity we assume non-degenerate Wi, ordered increasingly. The amplitudes 
Vn^ = {l\Nn) are the projections of the states |iV„) (generated by the operators 7V„ 
out of the vacuum) on the energy eigenstates {I \ of the scattering problem. 

In actual calculations one has to truncate the correlation matrix, since one can 
consider only a finite number r of operators. The generalized eigenvalue problem 

C(t)C^^)(i)=A«(t,io)C(to)C'''(t), fc=l,2,..,r. (4.7) 
allows one to find the energy eigenvalues Wi efficiently even for comparatively small 



values of t. In (4.7) we assume to < t, e.g. tg — 1. We require the 7V„ , n = 1, 2, .., r 
to be linearly independent, thus C(to) is regular and the generalized eigenvalue 
problem is well defined. Indeed, it may be transformed to a standard eigenvalue 
problem for the matrix 

C-Hto) C{t) C-Hto), (4.8) 

with identical eigenvalues A^'^^(t, to) and eigenvectors u^^\t) — {to)C'"''\t). The 
solution for the eigenvalues is 

X^''\t,to) = e-'^*-'"^^" . (4.9) 

Some approximations have to be done in the numerical calculations, among them 
the restriction to a finite set of lattice operators and therefore to a finite set of 
eigenvalues. However, as has been pointed out in [103], the correction due to the 
truncation is of 0(exp[— (t — tQ)Wr+i]) by a perturbation calculation. The range of 
t- values is dictated by statistical errors of the values of the correlation function and 
is usually between 1 and 8. 



4.2. Spectrum and Phase Shifts 

Consider quantum mechanics in a (d — l)-dimensional spatial box, with periodic 
(torus) boundary conditions with cubic symmetry. Assume that the interaction 
region is localized and smaller than the extension of the box. Solving the static 
Schrodinger equation outside the interaction region (V — 0, i.e. the Helmholtz 
equation) and expanding the solution into cubic {SO{d — 1, Z)) spherical harmonics 
one finds a relation between the discrete values of the momenta and the phase shift 
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at these points [97, 103, 99]. As a simple example we consider the scalar sector 
and neglect higher angular momenta (> 3). We find 

exp(2^^o(fc)) = + ' , (4.10) 

moo(K) - i 

3 kL 

moo{k) = Zoa{l]q'^)/qn^ where q= — , (4.11) 

ZTT 

i.e. TOoo(fc) = cot(5o(fc), the usual if-matrix. The function Zqo is a generalized 
zeta-function, and one has 

moo(^) = cot gTT for d=:2 , 

= tV E f°^-d=4. (4.12) 

This gives a relation 

(5(fc„) = - $(fc„) (4.13) 

between the quantized values of k and the phase shift. For d = 4 the functional form 
of ^{k) is determined numerically (cf. [149]). For d = 2 one has the simple form 
*i'(^ri) — knL/2 and thus 

2(5(fc„) + fc„L = 2?i7r (neZ), (4.14) 

and this quantization condition may be easily visualized for the cases without and 
with interaction (fig. |2^ ). Because of the periodic b.c. of the wave functions one 
has ^(0) — '^{L), ^'(0) = 'i'(L). Due to the (localized) interaction a non-zero phase 
shift affects the quantization of the momenta. 



Fig. 23. The quantization in d=l (a) without and (b) with interaction phase shift. 

These arguments and the equation ( [4.12| ) may be applied to a quantum field 
theory problem under the following restrictions and assumptions. 

— The interaction region is localized < L and, in particular, the size of the single 
particle states is smaller than the lattice: 1/m ^ L. This is due to that one 
wants a pair of two free particles as asymptotic states. 

— The polarization effects are controllable. Of course there will be interaction, 
e.g. self interaction of the light particles around the torus world. In fact this 
is perturbatively the leading contribution. This effect can be understood and 
accounted for. 
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— On is confined to the elastic regime 4 < s < 9 (for a 3-particle inelastic thresh- 
old, or 16 for a 4-particle threshold, depending on the quantum number of the 
channel) . 

— One understands the lattice artifacts (they are 0{a?)). This is because we work 
close to but never at the phase transition. These terms are due to the grained 
structure of the lattice and leading contributions may be studied already in the 
non-interacting gaussian model. 

For the determination of the energy spectrum one should consider correlation 
functions of a sufficiently large number of observables with the correct quantum 
numbers, capable to represent the eigenspace of scattering states [103]. 

To summarize, for a given set of couplings the determination of the phase shifts 
in the scattering sector (at total momentum zero) consists of the following steps. 

1. For given lattice size L determine the single particle state; check for polarization 
effects and get the mass m. 

2. Determine the correlation matrix C{t) for various values of t, solve the eigenvalue 
problem for this matrix, and obtain the spectrum from the exponential decay 



(4£). The dimension r of C cannot be chosen too large, because then the 
inversion becomes numerically unstable. We denote the energy eigenvalues by 

Wn{L). 

3. From the values Wn{L) one finds fc„ with help of the dispersion relation 



Wn = 2y/m^ + fc2. (4.15) 

This uses the fact that the 2-particle energy is the sum of the (back to back) 

outgoing single particle energies. 

4. From the values fc„ one obtains 5(fc„) from ( [1.13| ). 
For each lattice size L one obtains in this way several (typically the lowest 1-4) 
energy eigenvalues and thus pairs fc„, (5(fc„). The procedure is then repeated for 
other lattice sizes. Eventually one gets more and more values k and i5, filling the 
domain < fc < -\/5/4 (or \/3 if the 4-particle threshold is the first inelastic one). 
In fig. ^ the spectrum of free 2-particle states (vanishing phase shift) for the 4D 
system is shown. 



4.3. Models and Results 

Here I want to discuss two studies of scalar models and their results, concentrating 
on common features. 

4-3.1. 2D Resonance Model 

The authors [45, 46, 44] study a model, which describes two light particles </5 that 
couple to a heavier particle rj giving rise to resonating behaviour. For this one 
couples two Ising fields through a 3-point term, 

'S' = 5'ising(<^, Kv) + 5'ising(»/,K,,) + | ^ Vx^xi'^x-p. + Vx+fi) ■ (4.16) 

xeA,^^=l,2 

The values of the fields are restricted to { + 1,-1} and S'ising denotes the usual 
action of the Ising model. The 3-point term is introduced in a nonlocal manner. 
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Fig. 24. The free 2-particle energy spectrum (5 = 0) according to (4.13); the density increases 
with the volume (from [148]). 



because ip'^ = 1. The system is studied in the symmetric phase with regard to both 
nearest- neighbour couphng K^p and Kjj. 

The case g = corresponds to the situation of two independent Ising models, 
each with a 2nd order phase transition at Kc = ^ \n{l + ^/2) ~ 0.44068. In the scahng 
hmit the model describes an interacting boson with mass m = — log(tanh k) — 2k 
[77]. Reformulating the Ising model at the phase transition as a theory of non- 
interacting fermions, it was shown that the scattering matrix assumes the value — 1 
independent of the momentum [119]. 

In the coupled case (g > 0), if we identify rj and ip with particle states, the 
term proportional to g gives rise to transitions like rj — > ipip rendering rj a resonance 
in the (pip channel, when kinematically allowed. The phase diagram of this model 
has interesting structure with limiting cases including random external fields and 
random bond models. Its structure is discussed in more detail in [46]. 

The values of the couplings K^p, are adjusted such that — 0.19 and to^ = 
0.50. The system is studied for the case of two independent Ising models g = as 
well as for the values g=0.02, 0.04. The lattice sizes L x T considered were in the 
range L = 12 ... 60 and T = 60 . . . 100. Cluster updating was used, with typically 
200 K measurements per point. In the 2-particle channel 4-6 operators Nn were 
considered and 2-5 lowest eigenvalues determined from diagonalizing the correlation 
function at distances At = 2 ... 8. 

4.3.2. 0(4) $4 Model 

Zimmermann et al. [148, 149, 139] study the \ ^ oo limit of the 0(4) model, i.e. 
the nonlinear cr-model, with a small external field. As discussed, this is a theory 
with three degenerate light pions tt (would-be- Goldstone bosons) and a heavy a (the 
Higgs boson). 

The bare coupling parameters n and j were adjusted to have renormalized masses 
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771^ = 0.23 and ma — 0.70. The system was simulated with the cluster algorithm on 
lattices of size x T with L = 12 ... 32 and T = 32 ... 40 (40-50 K measurements per 
point). The authors considered 6 operators in the 2-particle channel and determined 
typically 2-5 lowest energy eigenvalues from the correlator at distances At = 2 ... 8. 

4-3.3. Single Particle States 

For the determination of the energy spectrum a precise knowledge of the single 
particle mass and related finite size effects is important. The operator of a 93 state 
with momentum 

P^^2ttv/L, v^-L/2 + 1,...,L/2 (4.17) 
may be represented through 

'P{Pu,t) = j-^^'fi{x,t)exp{ip,,x) . (4.18) 

X 

Its connected correlation function over temporal distance t decays exponentially 
cx exp (— |i|) defining the single particle energy i?^; in particular we have -Ejy=o = 

The observed mass, as compared to the "real" mass at vanishing lattice spacing 
and infinite volume, incorporates contributions from polarization due to self inter- 
action around the torus (decreasing exponentially with L) and lattice artifacts due 
to the finite ultra-violet cutoff (polynomial in the lattice constant a). 



Fig. 25. Results {Eu,Pv) as obtained from the analysis of the single particle propagator i n the 
2D Ising model [g = 0, L = 50, fr nm [4 ^]) The dashed curve represents the continuum d.r. (4.19) 
and the full curve the lattice d.r. (4.2C). 



Fig. |2^ is a plot of vs. p^, obtained for the pure Ising model {g = 0) and 
= 50. Comparing the values {E^,pjj) with the continuum spectral relation , 



(4.19) 
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one observes deviations 0((ap)^), as expected [103]. Using the dispersion relation 
for the lattice propagator of a gaussian particle with mass m, 

= arcosh(l — cospv + coshm) , (4.20) 

one finds excellent agreement with the measured values. Thus this particular single 
particle state shows little deviation from the free lattice form. These results for the 
2D Ising model holds for the situations g ^ Q as well and was confirmed in the 4D 
study [148, 149, 139], too. 

Let us, for the moment, stay with the 2D Ising model. It may be considered the 
limit of a $^ theory with infinite bare 4-point coupling A. We arc in the symmetric 
phase, thus we expect the leading interaction to be due to an effective 4-point term. 
The finiteness of the spatial volume allows self-interaction with one particle running 
around the torus. This leads to a finite volume correction to the particle mass pro- 
portional to \re~™'^ I^/Tj. In the 2D model excellent agreement with this behaviour 
was observed. This mass shift on finite lattices has also been confirmed in Monte 
Carlo simulations of the Ising model in 4 dimensions [108]. 

4.-3.4. Scattering Sector 

In the 2-particle channel one considers operators with total momentum zero and 
quantum numbers of the ?y, a set as complete as possible. These operators may be 
constructed from the p = single ry-operator, 

N^it) = r,{p = Q,t) = ^^Y.'^{x,t) (4.21) 

X 

as well as from the combination of two c^-operators with opposite momentum, 

Nn{t) = Y,fnipM-p,t)'fi{p,t) , (4.22) 
P 

<P{P, i) = E e''"^"-^^/ t) , (4.23) 

X 

where p and x have d — 1 (integer valued) components. The wave function in mo- 
mentum space fn{p) has been chosen in various ways. 

— Plane waves would be the cigenstates for the non-interacting case. This corre- 
sponds to choosing all integer values of jpj^. 



Up) = S{n,\p\) for d = 2, 
Up) = S{n, |p|2) for = 4 . 



2^,.. (4-24) 



However, one could also choose special forms like the periodic singular solutions 
of the Helmholtz equation 

/"(P) = (4-25) 

as has been tried in [149]. 
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The question here is which choice gives the best representation of the scattering 
eigenstates. Any "complete" set would do in principle, but in practice one is limited 
to a small number of operators and thus there may be "bad" and "better" choices. 
In the calculations discussed here no clear priority for one of the discussed forms 
was found. Both worked equally good; thus one favours the simpler form (4.24). 

These operators and products thereof are now measured for each configuration 
and the resulting correlation function (4.4) is computed for various values of t, 
diagonalized and the energy spectrum determined from the exponential decay (|4.9|). 



Fig. 26. The lowest energy levels for (a) the non-interacting gaussian model and (b) the Ising 
model; the decoupled ?7-state is also shown. The full lines denote curves due to phase shifts (a) 
5 = and (b) 5 = - f (fig. from [46]). 



Fig. 27. The phase shift obtained for the Ising model from the energy spectrum in fig. |26| (fig. 
from [46]). 

In fig. we compare the results for the energy levels for the 2D Ising with those 
for the free, non-interacting case. We find that the lowest level moves from the 
threshold (in the free case) to higher values in the Ising situation. In the figure we 
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also plot the results of the (decoupled) single rj state at to,, = 0.5. The level crossing 
therefore is fictitious and does not really occur, since the ipip channel is decoupled 
from the ry-channel {g = 0). We note that we have only a few energy eigenstates 
between the elastic threshold (at 0.38) and the inelastic threshold (at 0.76). 

The critical 2D Ising model has an S'-matrix equal to — 1 and we therefore expect 
a phase shift (Rising = — f (mod n) in the scaling regime. Fig. p7| shows the results 
for the phase shift S{k) determined from the data for the Ising model (fig. ^6|(b)) as 
a function of the dimensionless momentum k/rriip. 

Figs. ^ exhibit results of the 2D resonance model in the coupled situation, for 
g ^ 0- For most of the points the errors are smaller than the symbols, typically 0.5% 
of the energy value. The degeneracies in the energy spectrum have now disappeared, 
resulting in avoided level crossing with a gap that grows with g. 



Fig. 28. Left: The measured energy levels for the ipip channel for the coupling constant values 
g = 0.02, vs. the spatial lattice extension L. The full curves show the theoretical expectations for 
the values from the phase shift as discussed in the text. The dashed lines indicate the 2- and the 
4-particle thresholds. Right: The phase shift as determined from the energy levels (fig. from [46]). 



One should keep in mind that points at neighboured values of k may come from 
quite different values of L and different branches of energy levels. In order to empha- 
size this situation corresponding symbols denote the energy values on one branch, 
i.e. for one value of n. Each branch, from small to large L contributes to the whole 
range of k values; it starts at larger values of W (or k) with a phase shift close 
to 7r/2, passes through the resonance value at the plateau and approaches — 7r/2 
towards smaller k. The overall consistent behaviour is impressive. The high energy 
values, corresponding to large k may feel the lattice cutoff effects O(a^) as well as 
possible misrepresentation of the energy eigenstates by the considered operators. 
Within the accuracy of the data none of these effects seems to be a problem. 

The observed overall behaviour has a simple interpretation: A resonating phase 
shift superimposes the background Ising phase shift. The resonance parameters may 
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be determined from the phase shift hke from experimental data, e.g. by a i^-matrix 
fit or an effective range approximation. 



Fig. 29. The energy levels in both s-wave channels of the 4D 0(4) model for a particular choice 
of couplings (mir = 0.295, = 0.906); the dashed curves denote results due to renormalized 
perturbation theory (from [148]). 



Fig. 30. The isospin s-wave phase shift with the Higgs resonance determined from the energy 
spectrum in a non-perturbative Monte Carlo simulation (from [148]). 

Similar results were obtained in the study of the 4D 0(4) model. In 4 dimensions 
the numerical task is much more demanding and therefore one could not study as 
broad a range of lattice sizes as in the 2D case. However, the region of the first 
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avoided level crossing could be reached, as demonstrated in fig. Fig. ^ gives 
the isospin s-wave phase shift; this is the channel with the Higgs resonance. 

4.4. What Can We Learn? 

Comparison to Perturbation Theory 
The value of the partial width agrees with the Born term approximation using the 
value of the renormalized 3-point coupling [148, 139]. This implies that the situation 
is indeed described by an effective Lagrangian with tree-level expansion. However, 
the effective couplings have to be determined from the non-perturbative Monte Carlo 
calculation. 

4-4-^- Representation of Operators in the Scattering Channel 

In realistic situations we have to consider scattering between boundstates: In QCD 
the asymptotic states are bound states of the original quark fields, indeed a very 
demanding situation. One of the results of the purely bosonic scattering studies 
presented above concerns the quality of the representation of states through the 
lattice operators Nn{t). 

If the set of lattice operators Nn{t) is complete, we may relate the eigenvectors 
of the diagonalization (4.8) to the physical eigenvectors C^'^^t) oc v*^') [45]. Let us 



express the original approximation (4.5) to C{t) in terms of these eigenvectors. 



C„,„(t)=^C^')*C^)|v«|2e-*'^' , n,m=l,2...r. (4.26) 

In particular consider the diagonal element 

C^n»W=ElCi"l'|v«|^e-*^'. (4.27) 
1=1 

For a continuous energy spectrum, 

Cnn{t) = J dWpn{W)e-''^ , (4.28) 

and we identify |Cn ^ l^jv^'-' [^ c>c pn{Wi), the spectral density of the correlation function 
of operator Nn. We therefore find that |Cn^|^ indicates the relative weight of the 
contribution of operator 7V„ to the energy eigenstate 

Fig. P shows the relative weights |Cn^|^ for / = 1, 2, 3 and n = 1, 2, 3 for the 2D 
resonance model (case g=0.04). From the figure it can be seen, that the most im- 
portant contributions to an eigenstate come from lattice operators that have similar 
eigenenergies in their eigensystem. Compare the energy spectrum of fig.s |2^ and ^ 
with fig. |3^. For small L the lowest energy state Wi is close to the resonance energy. 
There the state is completely dominated by the ?7-operator A^i (circles in fig. p8| ), 
which describes an ?7-particle at rest. For increasing L, Wi decreases and approaches 
the 2-particle threshold. This manifests itself in a drastic decline of the contribu- 
tion of A^i , accompanied by an increase of the amplitudes for two (/3-particles at rest 
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Fig. 31. Contributions of the operators Ni,N2 and A^3 (circles, triangles, squares) to the energy 
eigenstates (a) |1) , (b) |2) and (c) |3}, displayed as a function of the lattice size L (fig. from [46]). 



(operator N2, triangles) and with relative unit-momentum (operator N3, squares). 
For lattice size larger than 20, the energy Wi is already too far below the resonance 
energy so that the energetically higher 77 at rest cannot contribute much to |1). 

In this way we may discuss the other levels as well. E.g. the contribution of the 
77-operator Ni (circles) to the second lowest state |2) assumes peak values where it 
crosses the resonance energy plateau. 

Also, with this consideration we can identify misrepresentation of energy eigen- 
states in terms of the considered operators. At higher energies one observes a shift 
of the weight factors |Cn''P towards higher operators Nn- If one has a situation, 
where the center of weight lies outside the scope of operators considered, more op- 
erators should be included for a good representation of that state [45] . The study 
of the representation of scattering states by the operators entering the correlation 
functions therefore sheds light on the problem of constructing good approximations 
for the energy eigenstates and it allows to control unwanted truncation effects. 



4- 4- 3- Lattice Artifacts 

In [44] the calculations for the Ising model were repeated for different coupling, 
corresponding to coarser lattice spacing. One found a stronger deviation from the 
expected phase shift at am^ = 0.5 than at amcp = 0.19, compatible with the ex- 
pected corrections O(a^). However, as mentioned earlier, the d.r. ( 4.15| ) gives the 



total energy of the asymptotic 2-particle state which (under the assumption of lo- 
calized interaction region) is just twice the energy of the outgoing particles. Now 
the O(a^) corrections of the single particle d.r. ha ve be en nicely describ ed by re- 



placing the continuum d.r. by the lattice relation ( 4.29 ). If one replaces ( 4.15 ) by 
the corresponding lattice expression, 

Wn — 2arcosh(l — cosfc„ -I- coshm) , (4.29) 

the data for Wn now produce slightly different values of A;„ and S, in much better 
agreement with a constant value of —Tr/2. One concludes, that the leading O(a^) 
corrections can be expressed by replacing the continuum d.r. by the lattice relation. 
This observation was confirmed in the study of the 4D model [148]. 
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In summary, it seems that Liischer's method to determine phase shifts from the 
volume dependence of the discrete energy spectrum works nicely, also for resonances. 
There are problems, where some of the restrictions cannot be met. In cases where 
more than one channel has non-vanishing phase shifts the determining equation in- 
volves matrices and the identification becomes problematic. Most likely one needs 
further continuity arguments like those used in the phase shift analyses of experi- 
ments. 

Conclusion 

During the last decade computer simulation have become an important tool in QFT. 
Within the realm of purely bosonic field theory efficient stochastic algorithms have 
been developed, that allow high precision results. Almost always, however, it is 
insufficient to just let the programs run and wait long enough until the answer 
stabilizes. Finite size and relaxation problems are persistent. In certain limits we 
may have an idea of the structure of an underlying (possibly: effective) field theory 
that will have unknown parameters and a limited domain of validity. In this way 
different effective models may describe (iiffer(;nt limiting cases of the microscopic 
theory. The nonperturbative Monte Carlo calculations can provide the unknown 
values of the parameters in ab initio calculations. This has made the results based 
on a finite volume theory (like finite size scaling from renormalization group, chiral 
perturbation theory, or the finite volume dependence in the situation of phase shifts, 
to mention the examples discussed here) the most reliable ones. 

The good agreement of Monte Carlo results with analytical expectations has 
opened the way for more sophisticated applications like using the Monte Carlo ap- 
proach as an "experimental" test laboratory for quantum field theory models. Ob- 
viously one cannot hope to prove something (maybe one can disprove some ideas) 
but one can hope to identify relatively quickly promising or most-likely-bad theories. 
An example for this approach is the study of Yukawa-models and chiral fermions 
(cf. the review in [27]). The more conservative reason for Monte Carlo studies, 
and still the main motivation, is the possibility to determine numbers in this direct, 
non-perturbative quantization scheme, that are up to now not at all obtainable with 
other methods. 
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